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Abstract 

We investigate positron annihilation in the electron g clS clS cl CcLSG study for many-body theory, 
in particular the optimized Fermi Hypernetted Chain (FHNC-EL) method. We examine several 
approximation schemes and show that one has to go up to the most sophisticated implementation 
of the theory available at the moment in order to get annihilation rates that agree reasonably well 
with experimental data. Even though there is basically just one number we look at, the electron- 
positron pair distribution function at zero distance, it is exactly this number that dictates how the 
full pair distribution behaves: In most cases, it falls off monotonously towards unity as the distance 
increases. Cases where the electron-positron pair distribution exhibits a dip are precursors to the 
formation of bound electron-positron pairs. The formation of electron-positron pairs is indicated 
by a divergence of the FHNC-EL equations; from this we can estimate the density regime where 
positrons must be localized. This occurs in our calculations in the range 9.4 ^ v s ^ 10, where r s 
is the dimensionless density parameter of the electron liquid. 

PACS numbers: 78.70.Bj,71.10.Ca 
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I. INTRODUCTION 



The process of electron-positron annihilation has been studied intensively for several 
decades. In recent years positron annihilation spectroscopy has been routinely used for 
studying the electronic structures of solids. As far as the two-body process is concerned, 
the appropriate theoretical framework is quantum-electrodynamics (QED). Differential cross 
sections and annihilation rates have been examined for two-particle systems like positronium 
in much detail, and can be found in standard textbooks Coincidence measurements 

of gamma emission give the angular correlation of annihilation radiation (ACAR), that 
yields information about the electron momenta. As the recent discovery of the electronically 

n 

stable bound positron-Li state showed [3(, positrons in contact with neutral atoms can also 
be an interesting few-body system. Since then the list of atoms binding positrons has 
become long; the underlying calculations are usually performed using either the stochastic 
variational method or the configuration integration method. 

Annihilation of positrons in matter adds a many-body aspect to the problem. The 
experimental annihilation rates are by now well established, we shall compare our results 
with the data measured by Weisberg and Berko 4] for alkali metals. In his pioneering 
work, Ferrell \^ gives intuitive formulas for the annihilation rates. The simplest "many- 
body" methods use single particle wave functions; in this case, the only true many-body 
effect is the Pauli exclusion principle acting between electrons, but many-body correlations 
between the interacting particles are not taken into account. Qualitatively, correlations 
can be introduced by applying enhancement factors as described by Brandt p. A popular 
method for electron structure calculations and potentially also for positron systems j3| is 
density functional theory (DFT) . In particular, there are numerous applications of DFT to 
defects in solids. In DFT one can write the annihilation rate in terms of the electron and 
positron densities, and an enhancement factor |8[ to account for the excess electron density 
near the positron, in other words, to describe electron-positron correlations. 

The enhancement factor has been the subject of many recent studies sun Once 
the enhancement factor is known, one can apply standard DFT under various approximate 
schemes, like the local density approximation (LDA), generalized gradient approximation 
(GGA) or the weighted density approximation (WDA). One can then also evaluate the 
partial annihilation rates due to valence and core electrons. The density functionals are 
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derivatives of known properties of electron gas or electron-positron mixtures, and their 
quality has been tested only in the case 

far there have not been many attempts to formulate a microscopic many-body theory that 
deals with an inhomogeneous electron gas. A first move to this direction was made by 
Stachowiak and Borohski ["L2I ]. who studied the case of a spherical inhomogeneity in jellium. 

Most of the many-body aspects of positron annihilation rates are reflected in a single 
number, namely the value of the electron-positron distribution function, g 1B (r), at the 
origin. (As a convention, we shall label all two-body quantities that involve one positron 
("impurity") and one electron ("background") with a superscript IB.) The annihilation rate 
of a positron in homogeneous electron gas can be written in the form [lHI ]. 

1 12 



t r 3 
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g 1B (0) x lO^ec" 1 , (1.1) 



where r s is the familiar dimensionless density parameter of the electron gas. The r s -factor in 
formula (jl.lj) is merely a geometric factor that takes into account the decreasing probability 
of finding an electron at the locations of the positron due to decreasing electron density. 
The "enhancement factor" g IB (0) accounts for electron-positron correlations. These can be 
strong in a metal, hence g 1B (0) can be large. The positron impurity is delocalized at low r s 
and cannot give rise to any appreciable local enhancement of the electron density. Instead, 
there is an increased probability of finding an electron near the positron: This tendency is 
visible only in the pair correlations, not in the density. This explains why annihilation rates 
computed using a homogeneous electron gas agree well with the experimental data up to 
r s ~ 5. On the other hand, the electron density enhancement around a localized positron 
can be included, for example, in the spirit of LDA in DFT, where the electron-positron 
distribution function computed for the homogeneous case is multiplied with the spatially 
varying electron density 14, lh, 3 • 

The calculation of g 1B (r) is a matter of many-body physics, and the problem of determin- 
ing g IB (0) is evidently an issue of short-ranged correlations. Quite appropriately, it was dealt 
with within perturbation theory by solving the electron-positron Bethe-Goldstone equation 
[sj]. Brown, Jackson, and Lowy 0, Q| have also examined short-ranged electron-electron 
correlations in a Bethe-Goldstone theory and pointed out the possibility of electron-positron 



pair formation at low electron densities (large r s ). The theory of positron annihilation has 
been developed further by Borohski, Szotek and Stachowiak jlj| and Rubaszek and Sta- 
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chowiak 



20 



21] , and it appears to be able to reproduce the observed annihilation rates in 



simple metals. 

With the advent of highly re-summed variational techniques, a new generation of calcu- 
lations containing vastly richer diagrammatic structures than Bethe-Goldstone calculations 
was possible. As a physically relevant paradigm for a fermion mixture and electron-hole 
liquids, positronic impurities and electron-positron mixtures have been studied quite ex- 

"quasi-boson" approx- 
imation for the wave function of the electronic background, which maps the formalism to 
an effective boson theory. Closest to our approach are the calculations by Lantto j^J and 
Saarela j^. Compared with the former one, the present work has an improved diagram- 
matic summation and in particular a more consistent treatment of the antisymmetry of the 
wave function. Lantto [24j also employs a simplified version of the Euler equation which 
corresponds to our FHNC//0 approximation to be discussed in section HV Al but the the 
energy and the structure functions are evaluated using the full FHNC equations and thus 
violates the identities (j2.14|) . ()2.24j) and (jH.fij) . Saarela [25] modified the Euler-Lagrange 
equation of the pair-distribution function of a charged Bose system by adding an ad hoc 
electron-electron potential, such that the equation reproduces the exact free-fermion dis- 
tribution in the limit of infinite density. Although oversimplified and phenomenological, 
this approach gives g IB (0) close to the values obtained by Stachowiak and Lach 26] and 



Borohski and Nieminen 



141 ] . Stachowiak et al. have used HNC theory in combination with 



a Hart ree-Fo ck-typ e approximation and a self-consistent perturbation of a Jastrow state 



a riartree-i H ocK-ty 



The present work should be considered as an exercise in basic microscopic many-body 
techniques. We shall first outline the most complete version of the optimized Fermi- 
hypernetted-chain theory (FHNC-EL) j^. We will pay special attention to the full func- 
tional optimization of correlation functions, which removes all ambiguity from the opti- 
mization process. The technical details of our theory for a one-component electron system 
are described in Ref. 0, a more recent application to 3 He which includes the optimization 
of triplet correlations and the calculation of proper elementary diagrams may be found in 



Ref. 



3J. We have also recently examined a completely analogous problem in helium liquids, 



namely the calculation of properties of 4 He impurities in 3 He |35j. We shall then discuss 
the impurity theory and derive the relevant Euler equations. In the following sections, we 
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will lead the reader through a sequence of plausible approximations in order to determine 
what it takes to get the physics right. We will show that the very simple approximation 
of a mixture of charged bosons gives a reasonably good agreement with experimental data. 
However, a bosonic theory is unsatisfactory per se, but it turns out that the first fermionic 
corrections make things worse and the agreement is lost. Finally, we will show that only the 
full fermion theory produces results which agree again with experimental data. 



II. OPTIMIZED FERMI HYPERNETTED CHAIN METHOD 



A. Variational Wave Function 



This section gives a brief survey of the variational theory of a bulk Fermi liquid; the 



reader is referred to Refs. 



33, 



34 for details of the theory, and the diagrammatic definition 



of all technical quantities. Since no confusion can arise for the time being, we will in this 
section not spell out the particle species. 

The Jastrow-Feenberg theory |36| for a Fermi liquid assumes a trial wave function of the 
form 



F(r x , . . . , r N ) 



F(n,...,r^)$ (l 
1 



N), 



exp 



u 2 (r i ,r j )+ u 3 (r i ,r j ,r k ) + ... 

l<i<j<N l<i<j<k<N 



(2.1) 
(2.2) 



$o(lj • • • > N) is a model wave function, normally a Slater-determinant of plane waves. The 
correlation functions M„(r l5 . . . , r n ) are made unique by imposing the cluster-property 



Wn(r x , 



r„) ^ as |r; 



oo . 



(2.3) 



The wave function ()2.2|) is not exact; one way to see this is by realizing that the nodes of 
the wave function ()2.2|) are identical to those of the model function $o(l ; • • • In the 

parlance of Monte Carlo simulations this would be called a fixed node approximation. 



B. Fermion HNC equations 

Two components are essential for the execution of the Jastrow-Feenberg variational the- 
ory: The first is the development of cluster-expansion and resummation methods for the 



5 



pair distribution function in the homogeneous case, 

N(N-l)fd 3 r 3 ...d 3 r N \V Q (l, 



9{r) 



|2 ' 



(2.4) 



p 2 f (Pr!... d?r N \* (l,---,N)\' 
where r = ryi = | ri — r2 1 ; spin-summations are tacitly implied. The second component 
of the theory is the optimization of the correlation functions by minimization of the total 
energy 



(*o|ff|*o) 



5 w n (ri, . . . , r„) (tf I *o> 







for the Hamiltonian 



H 



i=l 2 m l<i<j<N 



(2.5) 
(2.6) 



where, in our case v(r) = e 2 /r is the Coulomb interaction. 

Let us first turn to the pair distribution function g(r), specifically to the FHNC equations 
determining g(r) from a given pair correlation function U2{r). Pair correlations are the 
most important ones, and in the case of electrons one usually neglects triplet correlations 



altogether. Justification 
3D charged Bose gas 



or this stems from the study of triplet correlations in the 2D and 
37j and verified by Monte Carlo calculations j^, note that triplet 
correlations are not the same as propagator corrections [34] which have occasionally been 
confused with Feynman-Cohen backflow. 

The FHNC equations are a set of four configuration-space, and four momentum-space 
equations formulated in terms of "nodal" iVy(r) and "non-nodal" diagrams Xi$(r) that 
are, in turn, characterized by their exchange structure {ij} G {dd, de, ee, cc}. Input to 
the equations is the pair correlation function v&(r), the Slater exchange function £(x) = 
3 (sin x — xcos x)/x 3 , and a set of "elementary diagrams" E^(r) 3^] that must be calculated 
one by one. 

The coordinate-space equations are 



'dd(r) = X dd (r) + iVdd(r) 

= exp [u 2 (r) + N dd (r) + E dd (r)} - 1 , 

de(r) = [l + T dd (r)][N de (r) + E de (r)}-N de (r), 

X ee (r) = [1 + r dd (r)] \--L\r) + N cc (r) + E ee {r) 

+ [l + r dd (r))][iVd c (r)+£d c (r)] 2 , 

X cc (r) = -lv dd (r)L(r) + E cc (r) 



(2.7) 



N ee (r) 



where v is the degree of degeneracy of the single-particle states, L(r) = tikpr) — 
v[N cc (r) + E cc (r)], and &f is the Fermi wave number. The "nodal" quantities iVy(r) are 
constructed in momentum space according to 

X dd (k) 



N dd (k) 

N dc (k) 

N ee (k) 
N cc (k) 



1-X dc (k)\ - [l + X ec (A;)]x dd (£;) 
1 - X dc (k) - X dd (k) 



X dd (k) 



1-X dc (k)\ - [l + X cc (A;)]x dd (£;) 

X dd (£Q + 2X dc (£;)+X cc (£;)-l 
[l-X de (fc)] 2 -[l + X ee (fc)]X dd (A;) " 

-x cci k) m/v-x c M 



X de (k) 



X ee (k) 



(2.8) 



l-X cc (k) 

We have above used the convention of defining a dimensionless Fourier-transform as f(k) = 



[f(r)f(k)^pjd 3 rf(r) 



ikr 



The pair distribution function can then be constructed from the above quantities: 

g(r) = [1 + r dd (r)] + N ec (r) + E cc (r) + [1 + N de (r) + ,E; dc (r)] 2 | . (2.9) 

The static structure function is equally important than the pair distribution function and 
has the relatively simple form 



S(k) = 1 + [g(r) - if (k) 

l+X ee (k) 



[1 - X dc (k)} 2 - [1 + X cc (k)]X dd (k) 
For further reference, we also introduce the quantity 

1-X dc (k) 



S d (k) -- 
and note the relationship 

Tdd(k) [1 - X dc {k)f - [1 + X cc (A:)]X dd (A:) 
The representation ()2.9j) of g(r) contains an explicit factor 



[l-X dc (fc)] 2 -[l + X cc (A;)]X dd (fc) 
X dd (k) 



r dd (r) = exp [u 2 {r) + N dd {r) + E dd {r)} 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



This form is therefore the natural choice when working in coordinate space and focusing 
on the strong short-range correlation structure. On the other hand, when we consider the 
static structure function S(k), the expression ()2.10|) is the more useful one. 



The naive implementation |40|, |4l| of the FHNC equations, referred to as FHNC/ /0 ap- 
proximation, would suggest, in analogy to the boson theory, to start with the omission of 
the "elementary diagrams" E^(r) and include these, order by order, as quantitative improve- 
ments as the theory is moved to the next level. In the FHNC/ /n approximation one keeps 
"elementary" diagrams up to the n-point diagram. However [32], such a procedure violates 
the exact features 

X dc (k) = 0(k) as k -> 0+ , 

1 + X cc {k) = S F {k) + 0{k 2 ) as k^0 + . (2.14) 

These properties originate from the Pauli principle and are particularly important for the op- 
timization problem [42]. They imply the cancellation of "elementary" and "non-elementary" 
exchange diagrams; in other words there exist classes of so-called "elementary" exchange 
diagrams that must not be neglected. We have dealt with these diagrams in an approximate 
way, dubbed as "/ /C" -approximation, which has been described and justified in Ref. 



C. Background Energy calculation 



The first step in deriving the Euler equations is the calculation of the energy functional. 
An important manipulation is the use of the Jackson-Feenberg identity 



FV 2 F = i(V 2 F 2 + F 2 V 2 ) + ^F 2 [V, [V, InF]] - - 



V,F Z 



(2.15) 



which shows that the expectation value of the kinetic energy can be divided into three parts, 



(f) =T F - 



Nh 2 P 



5m 



d 3 r g(r)V 2 u 2 {r)+T w . 



(2.16) 



Here Tp is the kinetic energy of the free Fermi gas, and Tjp is a kinetic energy term that is 
solely due to exchanges. We can write this term as 



8m ^ 



h 2 ^(^olfVijv^F 2 ]]!^) _n 2 N 



8m 



J rf 3 rV 2 Pl (r) 



(2-17) 



($o|F 2 |$ ) 

With we mean (indicated by the subscript £) a gradient operator that acts on the 
Slater determinant only. Operationally, Eq. ()2.17j) is to be understood as follows: First 
one calculates — disregarding the fact that the density is uniform — a cluster expansion 
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of the one-body density pi(r). The operator Vf then differentiates the exchange lines only 
that are attached, in such a cluster expansion, to the external point. This replaces the 
incoming and outgoing exchange lines £(\r — /cp)£(|r — hp) at the reference point by 
(h 2 /8m) V 2 J(\r - n\ k F )£(\r - k F ). Finally, one takes the limit of the uniform system, 
and integrates over the configuration space of the last particle. 

Combining Eqs. (|2.16|) and (j2.17|) with the potential energy provides us with the starting 
point for further manipulations, 

l = f + U d3r 9(r)v JF (r) + ^, (2.18) 

where 

h 2 

Vj F (r) = v(r) - —V 2 u 2 {r) (2.19) 
Am 

is the Jackson-Feenberg effective interaction. 



D. Fermion Euler equations 

The formal manipulations to derive an Euler equation for the optimal pair correlations 
are almost identical to the ones carried out for bosons. The variation with respect to the 
pair correlation function consists of two terms: One comes from the variation with respect 
to the pair correlation function U2{r) appearing in the Jackson-Feenberg interaction fjp(r), 
and the second one is due to the variation of the pair distribution function with respect to 
u 2 (r) and the variation of T JF : 

The contribution from the first term on the right hand side of Eq. (J2.2(Jj) to g'(r) is calculated 
in complete analogy to the Bose case by replacing, in turn, each correlation line exp [ttzfaj)] — 
1 by exp [w2( r ij)] ^w{ r ij)- The second term is calculated recalling the graphical construction 
scheme of Tjp described above, and applying the same procedure to a graphical expansion of 
g(r). Thus, the contribution to g'{r) originating from Tjp is obtained by replacing in g(r), in 
turn, every connected pair of exchange lines l{rijkF)t{r ik kF) by (h 2 /8m)V 2 £( y r i jkF)£(r i kkF)- 
Following this construction scheme, one derives a set of eight linear equations, the FHNC- 
equations, corresponding to the eight FHNC equations ()2.7|h (|2.8J) j^ . in which the Jackson- 
Feenberg effective potential and the differentiated exchange functions act as driving terms. 
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To derive a form of the fermion Euler equations that is useful for a numerical implemen- 
tation, we write the Euler equation (|2.20j) in momentum space as 

1 



-t(k) [S(k)-1)+S'(k) = 0. 



(2.21) 



where t(k) = h 2 k 2 /2m. The S'{k) is a linear combination of the non-nodal quantities -XL (A;) 



S'(k) 



E 



dX-^k) 



(2.22) 



ij'e{dd,dc,cc} 

where the X»(k) are constructed from the non-nodal quantities X^(k) analogously to the 
construction of g'(r) from g(r) described above. Next, we define three effective interactions 
in the dd-, de- and ee-channels as 



V dd (k) = X' dd (k) - -t(k)X dd (k) : 

V dc (k) = X' de (k), 

V ee (k) = X' cc (k) + ^t(k)X cc (k). 



(2.23) 



The quantity V dd (k), defined in Eq. (|2.23jl . may be identified with the "direct interaction" of 



the Babu-Brown theory 



431 ] of the quasiparticle interaction. From Eqs. (|2.14|) . the effective 



interactions inherit the long- wavelength properties 

V de (k) =0(k) as 
V ee (k) = O(P) as 



k 
k 



o- 
o- 



(2.24) 



Using the representation (|2.10j) for the calculation of S'(k) via Eq. (|2.22j) and eliminating 
the X'^{k) in favor of the Vy(fc) lets us rewrite 

S'(k) = S 2 (k)V dd (k) + 2S(k)S d (k)V de (k) + S 2 d (k)V cc (k) + ~ t(k) [S 2 d (k) - S(k)] . (2.25) 

Inserting this expression for S'(k) in the Euler equation 1)2.21)1 lets us express S(k) in terms 
of the three effective interactions V^(k). 

The second step in the derivation of the Euler equation is to eliminate the pair correlation 
function U2{r) by using the FHNC equation (J2.7)) 



Am 



[l + r dd (r)]V 2 n 2 (r) 



^v 2 r dd r +- 

4m m 

h 2 



VJl + T dd (r) 



+ ^ [1 + T dd (r)] [V 2 N dd (r) + V 2 E dd (r) 



(2.26) 
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Using Eq. ([2.26)1 . one can rewrite the effective interactions Vy(r) in coordinate space entirely 
in terms of the distribution functions, the yet unspecified sets of elementary diagrams E^{r) 
and their "primed" counterparts E'^(r). The resulting equations are lengthy and not very 



33 and 



3J. For the purpose of comparison 



illuminating; they have been spelled out in Refs. 
with the impurity results, we display the coordinate space form of the direct interaction 

Ji 2 



V dd r) 



wi(r) 



[i + r dd (r)] 



m 
Am 



v(r) + ^-V 2 E dd (r) + E' dd [ 



VJl + r dd (r) 



+ T dd (r)wi(r) , 



V 2 iV dd (r) + iV dd (r) 



(2.27) 
(2.28) 



To calculate the effective interactions Vij(fc) ((ij) G {dd, de, ee, cc}) (or XUk)) and the 
fermion analog of the induced potential wj(r), one must also calculate the "primed" analogs 
of the "nodal" diagrams iVy(r). These quantities may be found in Refs. QandQ they are 
needed for the numerical optimization, but we will not need them in the further discussion. 
The dd-elementary diagrams have no special features; their omission can cause quantitative 
changes in the final answer, but does, unlike the de and ee "elementary" diagrams, not 
change the analytic structure and the properties of the solutions. 



III. IMPURITY CORRELATIONS 



In this section, and further on, we also must spell out the particle species, which may be an 
electron ( "background particle" , referred to by a superscript "B" ) or a positron ( "impurity 
particle", referred to by a superscript "I"). It is not necessary to label those quantities that 
were introduced in the last section and that refer only to background particles. We will keep 
the formulation general in the sense that the impurity mass mj is arbitrary, as well as the 
interaction between an "impurity" (positron) and a "background" (electron) particle. The 
Hamiltonian for the full system including the impurity is given by 

i2 N *2 N N 

H 1 = -£-V* - £ A-VJ + E -"(Ire - r,|) + E - r h \) ■ (3.1) 

1 i=i i=i j,k=l 

j<k 

As a convention, the impurity particle coordinate is i*o. In the present case of a positron in 
an electron gas, the impurity-background interaction v 1B (r) only differs in the overall sign 
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from the background-background interaction v(r), which is simply the repulsive Coulomb 
potential. The impurity mass mj is equal to the background mass m, but for a better insight 
into the problem it is helpful to keep mj. 

The formulation of the FHNC-EL equations for a single impurity follows essentially the 
same path as the formulation of the background equations, namely 

• Define a variational wave function 



• Derive a set of FHNC equations, 

• Derive the corresponding Euler equation using the "prime equation" technique, and 

• Reformulate the Euler equation in terms of distribution functions, thereby eliminating 
any reference to the correlation function u^ij-ij). 

Compared to finite-concentration mixtures, the derivation is simplified because there are no 
exchanges connected to impurity coordinates. 

A. FHNC equations for one impurity 

The (F)HNC technique is well established in earlier work, so there is no need to go 
through the details of the derivations here. The aspect that distinguishes fermions from 
bosons are the combinatorial rules and long- wavelength properties discussed above. In a 
mixture, exchanges can only take place between particles of the same species; this implies 
in the dilute limit that exchanges occur only between "background" particles. Moreover, in 
the dilute limit of a mixture only one impurity can occur in each diagrammatic quantity. 

Since the impurity cannot be involved in any exchange, we have only two FHNC equa- 
tions: The equations describing the parallel connections between external coordinates are 



M/ I B (0,l,...,iV)=exp -E4 B (ro,r,) Mh---,N), 



(3.2) 



L i=i 



*de(0 = [l + rE(r)l \E£(r)+N£(rj\-N£(r) 




(3.3) 



while the chain connections are best written in momentum space, 



N™(k) = (S(k) - S d (k)) X™(k) + (S d (k) - 1 - f dd (k)) X™(k) . 



(3.4) 



12 



From these quantities, we can construct the impurity-background distribution function 



9 lB (r) 



i + r£(r) 



1+Ng(r) + E!»(r 



dc ' 



(3.5) 



The long-wavelength properties corresponding to the identities (I2.14J) apply only for the 
background coordinates. Since the exchange structures of X d B (r) and X de (r) are the same, 
we have the long-wavelength limit 



X£(k) = O(k) as k -> + . 



(3.6) 



To abbreviate the equations, we found it convenient to define the quantity 

A few relations are useful in the derivation of the resulting equations: 
S lB (k) = [^{r)-l]*(k) = S{k)X%(k)+S d {k)X%W 

X l d B (k) 



= s d {k)x$(k) + ru{k)xg(k) = s d {k)x m (k) 



l + X ee (k) 



fdd(fc) 



S(k) l + X ee (k) 



(3.7) 



(3.8) 



B. Euler equations for one impurity 

For the determination of existence and stability of the solutions of the optimization 
problem, it is again useful to formulate the Euler equation in momentum space; the inclu- 
sion of the appropriate classes of "elementary" exchange diagrams always guarantees the 
proper short-distance behavior. Useful abbreviations are tj(k) = h 2 k 2 /2mi and the Feyn- 
man spectrum of the background, huAk) = t(k)/S(k). The formal Euler equation for the 
impurity-background correlations is |44|, in analogy to Eq. (|2.21j) . 



1 



^(tr^+tik^S^W + S'^ik) 



0. 



(3.9) 



The remaining manipulations are to carry out the "priming" operation on the impurity 
FHNC equations and to formulate the equations in a reasonably plausible form. This can 
be done in many ways, and the ultimate choice of the formulation depends to some extent 
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on the iteration path adopted for the numerical solution. Formally, we can define — in 
analogy to X m {k) introduced above — and to Eqs. (|2.23|) 

V IB (A;) = X' m {k) - l - (t(jfe) + t^X^k) (3.10) 

and rewrite the impurity Euler equation (|3.9J) as 

^<*> = - 2 itfJw (3 ' n) 

This representation of the Euler equation is formally identical to the Euler equation for 
impurities in Bose liquids. Of course, we still need to derive working formulas for calculating 
the quantity V m {k). 



C. Induced Interactions 



From the definitions ()3.10j) and (|3.7)1 . we can write formally 



v 1B (fc) = 0*0 + 



S(k) 



0*) + 



+ 



v£(k) + ^&v£(k). 

S{k) 1+X ee (k) 
t(k)S d (k)X™(k)X cc (k) 



(3.12) 



2 S(k) 1 + X ee (k) 

with (note that we deviate, for convenience and consistency with the definition (J3.1Uj) slightly 
from the definitions ()2.23j) ) 



V^(k) = Xg(k) - - (tr(k) + t(k)) X%(k) , 
Vg{k) = *2?(*) - ~ (ti(k) + t(k)) X™(k) . 



(3.13) 



The calculation of V™(r) is identical to the one for bosons and mixtures: 



i + r^(r 



+ 



dd 

2 t2 



v lB (r) + AV^ir 



IB, 



+ 



ft 



2m 2m i 



vvi + r&(r 



+ T d B d (r)<(r 



(3.14) 



with the induced interaction 



<(k) = N£(k) - - (tj(k) + t(k)) N£{k) 



(3.15) 
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and the elementary-diagram correction 

AV™(r) = 

Finally, we need 

= E'™(r) (l + r dd 
+ r™(r)N'™(r) + 



tf_ h 2 ^ 

Sm 8m/ 



IB, 



(3.16) 



h 2 Ti 2 \ 

^ + ^)V 2 X™(r) 
Smi 8m J 



(3.17) 



where the de-elementary diagrams must be chosen to guarantee the property V^(k) — > 
as k -> 0+, cf. Eq. (I3~fil . 

For the calculation of the remaining ingredients, it is most convenient to start with the 
effective interaction of correlated basis functions (CBF) theory 



v£(k) = r'£(k) - -(tj(k) + t(k))f™(k) . 



(3.18) 



Again, using background and impurity Euler equations, one can write this quantity as 



Vg(k) = - 



h(k) + 



t(k) 



1+X ee {k) 



r£(*) 



l + X ee (k) 



Vlf(k) + S IB (k)V dc (k) + f™(k)V cc (k) + -(*,(*) + t(k))X™(k) 



From this, we can obtain, for example, the induced interaction 

w?(k) = V™(k) - Vg(k) 



(3.19) 



(3.20) 



and, of course, r dd (k) which is needed for the rfe-equation, and N^(k) which is now easily 
obtained from 

= -\(ti(k) + t(k)) (s 1B (k) + fjg(*) + X™(k)) - Vg(k) - Vfe(k) . 
With this, we have derived a complete set of equations that can be solved by iteration. 



D. Impurity Energetics 



Another physical quantity of interest is the chemical potential of the positron impurity, 
which is the energy gained or lost by adding one impurity particle into the liquid, i.e. the 
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energy difference 



/i — E N+ i — E N 



(3.21) 



In the calculation of the impurity chemical potential from the definition (|3.21|) we must in- 
clude, besides the explicit terms containing impurity distribution functions, also the changes 
in the background distribution and correlation functions due to the presence of an impu- 
rity. These changes are of the order of 1 /N and therefore cause a change of order unity in 
the positron correlation energy. For brevity, we suppress here the contribution from triplet 
calculations; these corrections are already discussed in Ref. 
potential is, at the pair correlation level, given by 



35l The impurity chemical 



+ 



p I d 3 r 

P 2 n 



[g 



113/ 



l}v 



A < / d 6 r 



lB 'r)-g m {r) 



n 2 



h 2 h 2 \^, 
r - + 5~ 1 v u 



2„.IB/ 



(g(r) - l)v(r) - ^—g(r)V 2 u 2 {r) 
Am 



+ AT- 



.IF 



(3.22) 



where the A's in the second line of Eq. f!3.22j) indicate that we take the difference of the 
expressions calculated for the full system minus the same expression for the pure background. 

The pair correlations between background and impurity particles are determined by op- 
timization. Having a relationship between the pair distribution functions g l ^(r) and the pair 
correlation functions «2 ( r ) allows us to choose which one of these four quantities we consider 
to be the independent ones. The most convenient choice is to use the "dressed" correlation 
functions Tdd(^) and T d ^(r) since all other diagrammatic quantities can be defined in terms 
of these functions, whereas U2(r) appears explicitly only in the coordinate space equation 
()2.7)1 . In other words, we consider U2(r) as a functional of rdd(r), the impurity density p 1 , 
and p. As mentioned earlier, the consideration is further simplified by the fact that the 
impurity cannot be involved in any exchange. 

Since we need the chemical potential p) only to leading order in the impurity density, the 
optimization conditions for the background can be used to simplify the expression ()3.22|) for 
the chemical potential (again at the two-body correlation level): 



h 2 



h 2 

r2,) " ^( V i + V 2Mri,r 2 ) 



^A jy" rf 3 rirf 3 r 2 p(ri,r 2 

= ~\j d 3 r 1 d 3 r 2 p(r 1 ,r 2 )^(V 2 1 + V^)Au 2 (n, r 2 ) . 



(3.23) 



16 



In particular, since Tjp can be expressed entirely in terms of the T^{fij) and exchange 
functions, there is no rearrangement correction from this term. The change, Au 2 , is now 
expressed as a functional of the impurity quantities, 



Au 2 (r) = Au 2 \p I ,g m ,u m] 



d 3 r p 



j 5 u 2 {r 1 ,r 2 ) 
bp 1 



Furthermore using the HNC equation we find, again in HNC approximation, that 

j 5AT dd (n,r 2 ) 



Au 2 (ri,r 2 ) 



cfr p 1 



Sp 1 



(3.24) 



(3.25) 



The variations are carried out for fixed pair distribution functions. The change in the sum 
of nodal diagrams has a simple expression in the momentum space, 

!SN dd (k) 



P 



dp 1 



(*£(*))■ 



(3.26) 



The HNC equations ()2.7)1 and ()3.3|) are now used to eliminate the pair correlation func- 
tions u 2 {r) and u l 2 B (r) from the chemical potential ()3.22j) . Using the HNC equations and 
Eq. ()3.26|) one finds after some algebraic manipulations 



' p d 3 r 



f'HNC 



g m (r)v Lli (r) + 



n 2 



+ 



h 2 



2m/ 2m 



VvV B (r; 



+ 
+ 



d 3 k 



A J (2tt) 3 p 
1 
4 



[t(k)+tj(k)} S lB (k) N£(k) + E£{k) - N£(k) + E£(k) TZ(k) 



IB/ 



r IB, 



'-.IB, 



'-.IB, 



,/ " / '' t(k)[S(k)-l]\x™(kf 



(3.27) 



(2tt) 3 p 

The expression given above must be supplemented by corrections originating from 
"proper" elementary diagrams and triplet correlations (see Ref. |35J) if appropriate. 



IV. SIMPLIFICATIONS 



The FHNC-EL equations for the background, and even more for the impurity, are ad- 
mittedly complicated and of little appeal. The reason is that the exchange structure allows 
for many different ways of coupling. 

We recall, however, that the original motivation for deriving these equations is a symmet- 
ric treatment of short-and long- ranged correlations, and the stability criteria provided by 
the optimization. One is therefore tempted to reduce the FHNC-EL equations to a level that 
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contains the minimal amount of self-consistency and is, nevertheless, optimizeable. An ap- 
pealing feature of such a simplified theory is that the equations are hardly more complicated 
than those of the Bose theory; they are also readily applied to the much more demanding 
problem of inhomogeneous systems, where the solution of the full set of FHNC-EL equa- 



45| . In this connection, we need to 



tions is a rather unpleasant and still uncompleted task 
recall two things: First, positron annihilation rates are directly connected to short-ranged 
correlations, c.f. Eq. (jl.lj) . Second, it is far from straightforward to deduce short-ranged 
correlations in an inhomogeneous geometry from results on the homogeneous electron gas 

y. 



A. FHNC//0 for the background 

The FHNC//0 approximation is, by construction, the minimum approximation that is 
(a) correct for both short- and long-range correlations and (b) permits optimization. The 
FHNC/ /0 approximation for the background amounts to ignoring all de quantities and using 
1 + X cc {k) = S F (k), where 

{— — fc3 , k < 2kp 
1 k>2k F 

is the static structure function of the non-interacting Fermi gas. 

Consistent with this approximation, we leave out "elementary" diagrams; they may be 
put in at the end on a term-by-term basis. The FHNC equations (|2.7p . ()2.8|) collapse to 
three equations: 



A dd (r) = exp [u 2 (r) + N dd (r)} - 1 - A dd (r) , 

A dd (A;) = , (4.2) 

ddV ' 1 - S F (k)X dd (k) y } 

s(k) - Sp(k) 



1 - S F (k)X dd (k) 

Eqs. (J4.2j) are exact in the long-wavelength limit, but can be applied also at finite momenta. 
In the same approximation, we obtain from Eq. (J2.21j) 

S F (k) 



S(k) 
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V dd (r) 

Wi(k) 



[l + r dd (r)]v(r) + 



27.2 



h k 

4m 



m 
1 



S F (k) S(k) 



V^/l + r dd (r 

s F (k) + L 



+ T dd (r)wi(r) 



For further reference, we also mention that, in the same approximation, 



1 



S F (k) 



(4.3) 



(4.4) 



B. FHNC//O for the impurity 



Leaving the FHNC equations of section lill Al unchanged but identifying 1 + X cc (k) = 
S F (k) and S{k)/ S d {k) = S F (k) and ignoring V de (k) and V ee (k), the impurity Euler equations 
reduce to 



vim t(k) x™(k) 

S F (k) 2 S 2 F (k) 



[S F (k) - 1] 



(4.5) 



The induced interaction is 



w\ B (k) 



ti(k) + 



S F (k) 



S F (k) 



v£{k)--Mk)+t(k))xg{k) 



- V£(k). 



(4.6) 



The cfe'-equation does not simplify significantly, note especially that it is not legitimate to 
omit the "elementary" exchange diagrams. 

Let us finally also ignore all impurity de quantities. The Euler equation (|3.11|) remains 
the same, we only need to replace X m — > X™ and V IB — > V d l ® . We can write it also as 



S m (k) 



V d ^(r) = - 



m 



h{k) + huj(k) ' 
+ (l + T d B d (r)) v lB (r) + T™(r)wl B (r). 



(4.7) 
(4.8) 



The induced interaction is 



<(k) = -V£{k) 



S F (k) 



+ tj(k) 



mk) 



1 QlB (k] ( t{k) t(k) tj(k) tj(k) \ 

2 \S 2 (k) S 2 F (k) S(k) S F (k)J 



(4.9) 
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and finally 

?£(*> = 3- 

We can apply the above sequence of simplifications in two ways: One way is to use the 
FHNC//0 approximation for the background only, keeping all impurity equations intact. 
This might be a possible strategy for calculating annihilation rates in inhomogeneous ge- 
ometries, the good agreement between the g lB {0) values in this approximation with the 



471 ] provides encouragement that this is perhaps a pragmatic 



Monte Carlo data of Ortiz 
way to attack this problem. 

The second level of approximation is to also omit the cie-diagrams in the impurity equa- 



tions. We note t 



rat both strategies include the self-consistent summation of ring- and 



ladder- diagrams 4js,|4j|, both strategies also maintain "perfect screening" S IB (k) = 1. The 
FHNC//0 approximation provides for bulk electrons an accuracy that is only marginally 
worse than the full theory j3], and one might hope that the same is true for electron- 
positron mixtures. However, since the above approximation omits the self-consistent sum- 
mation of "elementary" exchange diagrams, the pair distribution function and the structure 
function become inconsistent. A simple approximation for g 1B (r) would be 

g 1B (r)^l + T™(r). (4.11) 

Adding the simplest exchange term would lead to a different approximation, 



^(r)« [l + rS(0J[l + tf(r)] , (4.12) 

where C(k) = T^(k) [Sp(k) — 1]. This is acceptable in a system with repulsive interactions 
since l + r^(r) is small at the origin. However, the opposite is true for the electron-positron 
case because both g 1B (r) and l+r™(r) become large at the origin, thus enhancing any small 
error that might be made in calculating the exchange corrections. We will see in the next 
section that this effect can be quite dramatic. 



V. ANNIHILATION RATES 

The central quantity of interest of our calculation is the positron annihilation rate (jl.lj) . 
hence, our primary question is what does it take to get r right ? To examine this question, 
we have carried out the following sequence of calculations of increasing complexity: 
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• Charged bosons: Positive impurity in a charged Bose gas. 

• FHNC/ /0 : The basic version of the FHNC-EL theory for both the background and 
the impurity. This means we use Eqs. ()4.2|) - ()4.4|) for the background, and Eqs. ()4.7j) - 
(j4.10j) for the impurity. The motivation for this approximation is that it is the simplest 
one that contains fermionic corrections, satisfies all exact short- and long-wavelength 
properties, and can therefore be optimized. 

• FHNC//0b: The full impurity FHNC-EL equations on the simplified background. The 
motivation for this approximation is that there is evidence that a simplified treatment 
of the background electrons could be appropriate, in fact, the energetics of the bulk 
electron gas predicted by this approximation is not significantly worse than the ener- 
getics produced by the full theory 4^ . However, electrons and positrons are strongly 
correlated and, thus, the same approximation might not be adequate for electron- 
positron correlations. Since the treatment of the impurity correlations is identical 
to the one of the full FHNC-EL theory to be described below, the comparison of 
the results from this calculation with the ones from the next more sophisticated one 
examines the sensitivity of the results to the background correlations. 

• FHNC/ /CO : The full solution of the FHNC-EL equations for the electrons and the 
positronic impurity, including corrections that guarantee the long-wavelength limit 
(J3.6j) but without proper elementary diagrams. 

• FHNC/ /C5 : As FHNC//C0, but including 4 th and 5 th order proper elementary dia- 
grams as described in Ref. 134 

Table |U lists a comparison of our results for the primary quantity, <? IB (0), in the simplest 
charged boson calculation (see also Ref. I25L our most sophisticated FHNC//C5 calcula- 
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24 



471 ] . At high densities, one 



tion, and a few important earlier calculations 
has reasonable agreement whereas the agreement becomes worse at lower densities. The 
comparison of different methods should serve to identify the source of these discrepancies. 



Among the various calculations, the variational Monte Carlo calculations of Ortiz |47[ have 
produced consistently the lowest values for g IB (0); incidentally, these values are quite close 
to our FHNC//0b results. 
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r s 


Kahana 


BSS 


Lantto 


SL 


Ortiz 


Bosons FHNC//C5 


1 






2.16 


2.29 


2.06 


2.40 


2.076 


2 


2.21 


3.76 


4.05 


3.96 


3.39 


4.55 


3.983 


3 


2.67 


6.91 


7.40 


7.29 


7.00 


8.25 


7.658 


4 


3.17 


13.6 


13.2 


13.6 


11.70 


14.71 


14.455 


5 






23.0 


24.2 


17.69 


26.08 


26.225 


6 








40.1 




46.18 


45.640 


8 








93.9 


64.00 


143.30 


126.063 



TABLE I: The table compares the primary quantity of interest, g IB (0), for the "charged boson" 
approximation (Bosons) and the FHNC//C5 approximation, with the values obtained by Kahana 
Borohski, Szotek and Stachowiak (BSS) ^(|, Lantto j^J, Stachowiak and Lach (SL) and 
Ortiz 



oron; 



Fig. ^ shows the annihilation rates r _1 computed using the approximations listed above. 
More appropriately, one should refer to these numbers as "partial annihilation rates" because 
we have neglected the effect of core electrons. Sob |5i| obtains, for example, 5% core 
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arrives at even 



contribution for Li, and 30 % for Cs, whereas more recent work 
larger core corrections. We conclude first that, in the experimentally accessible density 
range, the "charged boson" approximation is not as bad as one might expect. Inclusion of the 
simplest fermion corrections (FHNC/ /0) yields, on the other hand, far too high annihilation 
rates. This failure is rather surprising, because the approximation already includes the self- 
consistent summation of both ring and ladder diagrams 0,E]> and is thus beyond standard 
perturbative treatments that include only one of the classes of diagrams or an incomplete 
self-consistency scheme. 

The main reason for the failure of the FHNC//0 approximation can be traced to the way 
exchange corrections were included (see the discussion in the end of section HV]) . On the level 
FHNC/ /0b this problem is solved by including the full FHNC summation for the impurity 
correlations, while leaving the background correlations at the FHNC/ /0 level. Indeed, there 
is some improvement. The annihilation rates come out lower than those predicted by the 
more sophisticated approaches to be discussed below. But, in fact, when core corrections 
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FIG. 1: Experimental and theoretical annihilation rates for various materials. T 



re experimental 
I54L without 



points are, from left to right, for the metals Al, Zn, Mg, Li, Na, K, Rb, and Cs 
core corrections. The theoretical curves correspond to various levels of implementation of the 
(F)HNC-EL theory as labeled in the figure. In the level of sophistication the order from simple to 
complex is: Charged bosons, FHNC//0, FHNC//0b, FHNC//C0, FHNC//C5 (see main text for 
definitions). The positronium decay rate is (500 ps) _1 = 2 x 10 9 s . 

are included, agreement of the FHNC//0b approximation with experiments is not too bad. 
The full FHNC summation FHNC/ /CO brings us back up, and ultimately the inclusion of 
"proper" elementary diagrams on level FHNC/ /C5 produces reasonable agreement with the 
experimental annihilation rates when core corrections are omitted. 

Evidently g 1B (0) is extremely sensitive to the description of the background electron 
liquid, as can be seen in the large difference between the full FHNC/ /0, the FHNC/ /Ob and 
the FHNC/ /CO results. We conclude therefore that any simplistic treatment of background 
correlations should be viewed with much reservation. 

In discussing the possible theoretical refinements one should finally comment on the 
importance of triplet correlations. Since the proper elementary diagrams have an impact, 
one might expect the same to be true for triplets as well. In the case of helium liquids, 
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-0.40 



-0.45 - 



P? -0.50 I- 



— 



-0.55 



-0.60 



-0.65 




FIG. 2: The calculated positron correlation energies as function of r s in the FHNC//C0 (dashed 
line) and the FHNC//C5 (solid line) approximation. Also shown are the theoretical results of 
24], Arponen and Pajanne (AP) s|, and Borohski and Stachowiak (BS) jsij ] - 



Lantto 



44. 



triplets have been mostly described in the "convolution approximation" outlined in Ref. 
Recently, the treatment has also been found to be adequate for a problem analogous to the 
present one, namely the calculation of the energetics of a 4 He impurity in 3 He 13511 . The 
next systematic improvement to the method of Ref. Q has also been examined js^. but 
for helium liquids it was found to be insignificant. Unfortunately, the same is not true for 
positronic impurities in the electron gas. We have tried the same approach here and found 
that the "convolution approximation" is manifestly inadequate for this case. It appears that 
the derivation of triplet HNC equations [56] is necessary for a satisfactory treatment. We are 
not aware of a successful implementation of this program for bosons, much less for fermions 
or mixtures, and must leave this question open for the time being. The variational Monte 
Carlo calculations of Ortiz [47] support our view that triplet correlations are unimportant. 

From now on we show only the FHNC/ /C5 result, and occasionally the FHNC/ /0 result 
for reference. The next quantity of interest is the positron correlation energy (or chemical 
potential) (|H.21J) . Our results are shown in Fig. they are consistent with earlier calcula- 
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tions of Arponen and Pajanne [9| and Lantto 



9] 



24|, although our results are somewhat higher. 



The present results also show a drop of the chemical potential towards lower densities; a 
similar drop has been reported by Hodges and Stott 57]. 

It has been argued that, in the low-density limit (large r s ), the correlation energy should 
go towards the binding energy of the free positronium of -0.524 eV. Evidently, our calculation 
does not show such an asymptotic behavior and we hasten to point out that it should not: 
The Jastrow-Feenberg function used in Eqs. (|2.2j) and (|3.2j) describes a state in which the 
positron is delocalized, so it is not a valid description of one positronium atom and N free 
electrons. Hence, the present theory is inapplicable to a system with bound electron-positron 
pairs. It has, however, the desirable feature that the Euler equation ceases to have a solution 
at the point where system cannot be in the state described by the trial wave function. In 
the present case, this happens at the density where the positron must be localized. This 
feature of the theory does not preclude that the free positronium is the physical low-density 
limit; the theory simply makes no statement about this limit. 

The appearance of a positronium atom is, on one hand, signaled by a very large g lB (0). 
Moreover, recall that we solve the Euler equation with the boundary condition g m {r) — > 1, 



as r — ► oo. If there were a bound state, then \Jg IB {r) should have a node at some finite 
distance. Of course — as argued above — the theory will not allow us to reach this point, 
but one should expect precursor phenomena as one approaches the density of positronium 
formation. Fig. |3]shows indeed both effects: As we decrease the density, the g m {0) increases. 
In addition to this, a "dip" forms at an intermediate distance. At the highest r s value shown 
in Fig. EJ r s = 9, we have g 1B (r) as low as 0.042. With decreasing density the dip drops 
rapidly, and finally the pair distribution function becomes unphysically negative at r s « 9.4. 
Also g 1B (0) diverges, at r s = 10, but because of both diagrammatic approximations and 
numerical limitations one cannot expect these instabilities to show up at exactly the same 
density. We conclude that our calculations predict positronium formation at 9.4 < r s < 10. 
For comparison, in the mid 70 's Lowy and Jackso n |18l| found the positronium limit at 
r s > 6.2, while a decade later Pietilainen and Kallio |23J] obtained r s > 8. 

Slightly different information is contained in the static electron-positron structure func- 
tion S IB (k), depicted in Fig. In the long wavelength limit, charge conservation implies 
S' IB (0+) = 1; this property is a rigorous feature of the optimized structure function in any 
level of the FHNC-EL approximations. As the density decreases, we see that S lB (k) de- 
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1000.000 




r (a r s ) 

FIG. 3: The figure shows our calculated electron-positron pair distribution function g lB (r) for 
r s = 1, 2, . . . , 9. The function with the lowest value of <? IB (0) corresponds to the highest electron 
density. Notice the logarithmic scale in g lB (r). 

velops a peak at k ~ 3.2(a r s ) _1 . Such a peak reflects an oscillatory structure in the pair 
distribution function. 

VI. SUMMARY 

We have presented in this work calculations of electron-positron correlations in simple 
metals. Technically, we have used the most complete summation of diagrams within the 
Jastrow-Feenberg theory. We have shown that such a highly summed theory is indeed nec- 
essary for a reasonably reliable prediction of the relevant quantities. The omission of triplet 
correlations is an unsatisfactory point, but progress towards a fully consistent summation of 
three-body (F)HNC equations and the corresponding Euler equation is not in sight. In any 
event, we have demonstrated that one might have to do even more, but one must definitely 
not do less than was done in our work to obtain a conclusive answer. 

Our results are, from a pragmatic point of view, somewhat disappointing in the sense 
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2 4 6 8 10 

k ((aor s )-i) 

FIG. 4: The figure shows our calculated electron-positron structure function 5 IB (/c) for r s = 
2,3, ... ,9. The lowest curve corresponds to the lowest r s (highest electron density). 

that the system is evidently much more complicated than the bulk electron gas at metallic 
densities. Bulk electrons at metallic densities are one of the simplest systems to be treated 
within microscopic many-body theory, and simple approximation provide already reasonable 
agreement with exact results. The reason for this problem is evidently the large value of 
the pair-distribution function at the origin, which causes poor convergence of diagrammatic 
summations. On the other hand, the FHNC/ /Ob approximation seems to provide reasonable 
agreement with experiments when core-corrections are included. With sufficient caution, 
since better agreement with experiments does not indicate a better theoretical treatment, 
one might be able to use this version of the theory in a non-uniform environment and obtain 
reasonable annihilation rates for positrons in simple metal surfaces j^]. 
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